incorporates the normal state ARPES dispersion, and non- 
magnetic impurity scattering in the presence of a d x2y2 or- 
der parameter. We have shown, that it is extremely im- 
portant to take into account the finite resolution of the an- 
alyzer in both momentum and frequency to obtain real- 
istic results. In particular, the broadening by finite mo- 
mentum resolution leads to a finite region of 'gapless- 
ness' around the true node of the d 2 2 OP, if the shift of 

x —y 

the leading edge of the spectrum, A sh , with respect to a 
non-superconducting reference sample (in electrical con- 
tact with the superconductor) is taken as the measure of the 
gap. This shift vanishes faster than linear in the Fermi sur- 
face angle (]) around the d r2 _ v2 node region, even though the 

x y 

true OP itself goes to zero linearly. 

Non-magnetic impurities amplify the effect of the fi- 
nite resolution. Additional spectral weight around co = 
is introduced, leading to a further reduction of the shift 
A sh . At sufficiently large n ( . « 0.02 — 0.05, this results in 
an extended region of apparent 'gaplessness' of 8(|) = ±7 
around the true node for a large range of moderate to strong 
impurity potential strengths. The zero frequency reso- 
nance responsible for a large residual DOS iV res at £ F man- 
ifests itself in the ARPES spectra in terms of a slope prob- 
lem (particularly prominent in the large gap region around 
[71,0]), i.e., the slope of the ARPES spectra and that of the 
non-superconducting reference sample are far from par- 
allel when the intensities are suitably normalized. There 
is a low frequency tail in the spectra reflecting the spec- 
tral weight from the impurity scattering near the resonance 
limit. We also demonstrated that this spectral weight could 
be identified by looking at the effect of particle-hole mix- 
ing for k values which are outside the normal state Fermi 
surface. In this case, the incoherent spectral weight around 
co = leads to a second peak in the ARPES spectra with a 
weight comparable to that of the quasiparticle part. 

One point about the impurity related parameters used 
in our calculations has to be born in mind: for n. > 0.02, 



|1/m| < 0.8, our model also leads to a significant reduc- 
tion in T (see Ref. 4) up to 50% at n. = 0.05. Hence, if 
an extended region of nodes is observed in an ARPES ex- 
periment, this should be accompanied by a substantial T 
reduction compared to a pure sample where such a 'node- 
less' region is absent. If this is not the case, it might well be 
that the sample surfaces have a higher impurity concentra- 
tion than the bulk, so that the T suppression is determined 
by a different n. than the one that affects the ARPES spec- 
tra. Note also, that NMR experiments on Bi2212 crystals 
grown by the same group as the ones used in the ARPES 
experiments Ref. 3, have indicated the presence of a large 
N without explicitly doped impurities. 20 This result is 
somewhat puzzling, since the reported value of N (24% 
of the normal DOS) would imply a T reduction of at least 
20% (according to our model calculations), but the T of 
the sample was still 86K, which is within 10% of the best 
Bi2212 values in the literature. Nevertheless, if this resid- 
ual DOS is due to non-magnetic impurity scattering, one 
would expect the presence of an extended nodeless region 
in ARPES measurements on those samples. 

A final word about resolution effects: Improved fre- 
quency resolution would be beneficial for the understand- 
ing of the ARPES line shape around (71,0), since due to 
the proximity of the saddle point, momentum broadening 
is very weak in this region, and one might be able to ob- 
serve more features of the true spectral function. On the 
other hand, better momentum resolution would definitely 
help to improve the accuracy in determining the size of the 
gap since in the most interesting region around the d x2 _ y2 
node, momentum resolution broadening is clearly the lim- 
iting factor due the rapid quasiparticle dispersion. 
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if the underlying spectral function deviates substantially 
from a single-peak structure. 
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Fig. 7. The superconducting gap on the Fermi surface as a 
function of the angle (ft. The true OP A is compared with A sh , 
calculated for l/u = -0.4 and different impurity concentra- 
tions n. 

uated for l/u = 0.4, n = 0.05. 



0.0,0.02,0.05. The thin solid line shows A sh eval- 



In Fig. 7, we show the angle dependence of A sh on the 
Fermi surface for l/u = —0.4 . As compared to the true 
OP, the shift A sh extracted from the EDC clearly underes- 
timates the value of A already in the pure case, an effect 
which is strongly enhanced by the broadening caused by 
impurity scattering. Furthermore, the impurities lead to a 
'gapless' region around § =45 of size up to 8(|) = ±7, where 
the shift A sh is essentially zero. Such behavior has been 
reported earlier for some samples measured by ARPES, 

1 Q 

while it is not seen in others. Our results strongly suggest 
that it is the sample purity which is responsible for these 
differences. A non-linear dependence on (]) of the OP itself 
is not necessary to explain the flatness around (]) = 45. 

The thin line in Fig. 7, was obtained for l/u = 0A,n. = 
0.05, and shows (in comparison to the curve for l/u = 
-0.4) that the shift A sh in the large gap region can be 
strongly dependent on the size and sign of l/u, while the 
curves essentially overlap in the small gap region around 
(]) = 45. This effect is related to the severity of the slope 
problem shown in Fig. 6, and can be traced back to the vari- 
ation of the peak position in the self-energies displayed in 
Fig. 4 as a function of l/u. 

Finally, we would like to illustrate another possible 
way to identify strong impurity scattering in ARPES ex- 
periments. Due to particle-hole mixing, the BCS spectral 
function Eq. (2) predicts that the ARPES intensity should 
not vanish abruptly when the normal state Fermi surface 
is crossed, but rather, a gradual decrease should occur, ac- 
cording to the dependence v k = (1 -t, k /E k )/2of the weight 
in the electron part of A(k,co). At the same time, the peak 
in the spectral function approaches the Fermi energy only 
up to the minimum distance at E, = —A, , and then dis- 

K F K F 

perses back to larger binding energy as k crosses the Fermi 
surface. Precisely this feature of BCS theory has recently 
been observed for the first time in Bi2212 19 making the 
case for the spectral function interpretation of the ARPES 
experiments on Bi2212 even stronger. 
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Fig. 8. Particle-hole mixing: Evolution of the EDC evaluated 
at 1 /m = 0. 1 , n . = 0.05 (close to resonance) as the quasiparticle 
energy moves through the Fermi energy, Z, k = 0.0,0.1,0.2,0.3. 
The lines indicate the dispersion of the electronic quasiparti- 
cle away from E F as the Fermi energy is crossed. The arrows 
mark the incoherent, impurity induced peak near CO = which 
becomes more distinct as £, k grows. 



Here, we demonstrate that impurity scattering may 
leave another unique fingerprint in the ARPES spectra 
when investigating the effect of the particle-hole mixing. 
The point is that as the Fermi surface is crossed, and the 
strength of the 'quasiparticle intensity', roughly given by 
v R = (1 —E, k /E k )/2, weakens rapidly, the spectral weight 
introduced by the impurity self-energy in the vicinity of 
co = decreases much slower. Hence, the incoherent part 
can show up as a second peak as soon as the coherent in- 
tensity (quasiparticle peak) is sufficiently small. This phe- 
nomena is demonstrated in Fig. 8 for a scattering strength 
l/u = 0.1 (near resonance, where the effect is strongest), 
and impurity concentration n. = 0.05. One clearly ob- 
serves how the quasiparticle peak disperses away from e F 
as the Fermi surface is crossed, while the low-frequency 
tail present for ^ k = progressively turns into a separate 
peak as cj k increases. We have plotted the EDCs without a 
normalization to show the drop in signal strength. Never- 
theless, the intensity might still be enough to be observed. 

V. CONCLUSION 

We have presented an analysis of ARPES data in the su- 
perconducting state of the high T cuprates based on a sim- 
ple phenomenological weak-coupling BCS model, which 
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Fig. 5. Spectral function, and the corresponding EDC calcu- 
lated using (1) for l/u = -0.4 and n ; = 0.02,0.05 at four dif- 
ferent angles on the Fermi surface, (]> — 0,35,40,45 (from top 
to bottom). First column: the unaveraged spectral functions. 
Second column: EDCs for n. — 0.02. Third column: EDCsfor 
n i — 0.02. The dotted line in columns two and three represents 
the resolution broadened Fermi function N R (u>, T). 

At (]) = 0, where the OP reaches its maximum, the shift 
is also strongly reduced due the progressive broadening of 
the EDC by the impurity scattering, but at n. = 0.05, a new 
effect appears: the leading edge slope of the EDC does not 
match the slope of the resolution function N R ((t>, T). This 
behavior is a consequence of the enhanced spectral weight 
around zero frequency which also induces the large N . 
The slope problem is most critical for scattering strengths 
around the resonance limit. 

Fig. 6 addresses this issue in more detail. Here, we plot 
the EDCs at (]) = 0, n ( . = 0.0, 0.05 for various values of 1 /u. 
For rather weak scattering (l/u = 1.0,-1.0), the slopes 
of the EDCs and the resolution function are pretty much 
parallel. However, in the strong scattering limit, mostly 
pronounced at l/u = 0.0, a low frequency tail in the EDC 
makes the peak substantially asymmetric, and results in the 
above mentioned slope problem. This effect could be a 
way to identify strong impurity scattering in ARPES ex- 
periments. We emphasize however, that it would be seen 
most easily in regions where the OP is large, i.e., around 

(71,0). 

In fact, a similar behavior has recently been reported 
from ARPES measurements on electron irradiated samples 
of Bi2212. 16 The electron irradiation is believed to create 



Frenkel-type point defects by displacing oxygen atoms in 
the Cu0 2 plane. This type of defect apparently does not 

i n 

lead to a Curie tail in the magnetic susceptibility, hence 
it seems to represent a non-magnetic potential scattering 
center. The spectrum shown in Ref. 16 (measured in the 
large gap region) indeed seems to have acquired a low- 
frequency tail as expected from the theory. A more de- 
tailed study would be desirable to unambiguously identify 
the slope problem. 
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Fig. 6. EDCs for various values of 1 /u and n ; . = 0.05 at the gap 
maximum (]) = 0, compared to the pure spectra, and the reso- 
lution broadened Fermi function N R ((0, T) . The slope problem 
is obvious close to resonance (l/u = 0.0, —0.4). 

Fig. 6 also clearly illustrates the problems one encoun- 
ters when trying to extract the OP from the leading edge 
shift A sh . Even though the OP is precisely the same for all 
spectra in Fig. 6, the value one extracts from A sh have a 
large variation. Such a slope problem will always occur 
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negative to positive frequencies as 1 ju is varied from 0.4 
to -0.4. This structure leads to the small peak in A(k F ,(0) 
provided it is away from z p . When it is at £ F , the case 
which corresponds to resonant scattering, the self-energy 
peak coincides with the quasiparticle peak positioned at cj k 
leading to a single structure in A(k p ,(a). 




-0.4 -0.2 0.0 0.2 0.4 

CO 

Fig. 3. Dependence of A(k p ,§) on the strength of the impurity 
potential u for n j — 0.05. Spectra are shown in the normal as 
well as in the superconducting state (at the node (]) = 45, and 
the maximum gap = 0). 

At (]) = (the large gap region) the spectra consists 
of the two peaks positioned at ±A k as expected from the 
BCS form of A(k,co). The impurity scattering substantially 
broadens the peaks, and it also leads to spectral weight at 
zero frequency, which is responsible for the finite Af. es . The 
weight around co = depends again on the value of 1 ju, for 



positive (negative) 1 ju it is concentrated at positive (nega- 
tive) CO. 

One further point worth mentioning is the fact that the 
impurity scattering violates the approximate sum rule stat- 
ing that n(k F ) should be independent of temperature above 
and below T ? It is evident from Fig. 3 that compared to 
its value in the normal state (= 0.5 by definition), n(k p ) in 
the superconducting can grow or decrease depending on 
the value of 1 ju. The effect is particularly severe at (]) = 45. 
Here, of course this is due to the strong particle-hole asym- 
metry introduced by the impurity scattering and the normal 
state dispersion, under which condition the sum rule is not 
applicable. 9 It might be that this effect could shed some 
light on the puzzling results from Ref. 15, where a large 
change in the ARPES spectral weight was observed upon 
cooling through T . 
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Fig. 4. Dependence of the imaginary part of the self-energy 
component Im{E n } on the strength of the impurity potential 
wfor n. = 0.02,0.05. 

Next, we look at the effect of the impurity scattering 
on the EDCs as calculated using (1). Fig. 5 illustrates the 
influence of the scattering on the shift A sh , at lju = —0.4 
for four angles (]) = 0,35,40,45 on the Fermi surface, and 
at two values of n j = 0.02,0.05. For comparison, the raw 
spectral functions are also shown in column one. At the 
node, (]) = 45 the impurities have virtually no influence on 
the EDC as evident by comparing with the spectrum in Fig. 
1 . However, at cf) = 40, increasing the impurity concentra- 
tion leads to a progressive closening of the gap measured 
in terms of the shift: at n. = 0.05, A , = 0. The same effect 

i ' sh 

is seen further away from the node at cb = 35. Here, A sh is 
still finite at n. = 0.05, but is substantially reduced from its 
value at n. = 0.02. 
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large gap region, the § dependence of A sh pretty much fol- 
lows that of the true OP, however, around the node, there 
is a big discrepancy coming from the finite k resolution: 
Whereas the d x2 _ y2 OP goes to zero linearly, A sh takes a 
super-linear dependence, and becomes zero already at an 
angle of (]) « 42, producing an apparent extended region of 
nodes around (]) = 45. Also, the gap is underestimated by 
a factor of at least 1.4 even in the large gap region which 
is of course obvious, since the true gap is given not by the 
shift but by the position of the peak in the EDC. This ex- 
ample illustrates that the k dependence of the gap around 
the node region is very difficult to quantify using the shift 
method already in a simple BCS model. In the following 
section, we shall show that impurity effects can even fur- 
ther blur the true value of the gap through lifetime effects. 

IV. IMPURITY EFFECTS 

As mentioned in the introduction, non-magnetic impuri- 
ties have a severe pair-breaking effect in a d wave super- 
conductor. In this section, we shall investigate their ef- 
fect on the spectral function, the EDC, and the values A sh . 
We use a short-range potential V.(r) = «8(r — r ; ), and ap- 
ply the self-consistent t-matrix approximation. 7 This ap- 
proach has recently been shown to yield accurate results in 
the dilute impurity limit. 14 Working in particle-hole space, 
the t-matrix T satisfies the following Lippmann-Schwinger 
equation (quantities with a hat represent matrices) 



r(co) = u 



g 3 + g 3 r(co)- J2^(k,co) 



(3) 



Here we introduced the propagator 



?(k,co) 



coa + A k a,+^ k o 3 



co z 



■ A 2 - E 2 

^k Sk 



(4) 



where co = co-E Q ,A k 



particle energy. The self-energy (like all other matrices) 
is expanded in terms of the identity and Pauli matrices 
G • • • G 3 as £ = ^IjGj, and is given by E = nJ, n ; the impu- 
rity concentration, and A the volume. For the d x2 _ y 2 OP, the 
off-diagonal self-energy vanishes by symmetry, and from 
(3), the non-zero components are 



n.G„ 
i o 



nj±ju -G 3 ) 



-0 (1/M _ G3 )2_ G 2' 3 (1 /m _G 3 )2- 



(5) 



where G^co) = A~ L k g ( .(k,co). As outlined previously, 
we take into account the full structure of the t-matrix, 
i.e., we do not assume particle-hole symmetry. The self- 
consistency equations (5) are solved numerically on a d x d 
k-space grid with d < 4000. The spectral function is then 
evaluated from 



A(k,C0) 



sgnco 

K 



Im{g n (k,C0)} : 



(6) 



and the EDC is obtained using Eq. (1). 

As shown earlier, 4 the effect of the non-trivial DOS is 
rather dramatic, in particular if a vHs is close to e F , as in 
the present case: (i) The s-wave scattering phase shift 8 
acquires a strong frequency dependence already in the nor- 
mal state, which also reflects itself in the scattering rate, 
hence resistivity. Furthermore, the dependence on u is 
highly non-trivial, (ii) In the self-energy Eqs. (5), the 
cotangent c = cot8 Q (being frequency independent in case 
of a constant normal DOS) that is usually used to parame- 
terize the scattering strength, is replaced by the (now fre- 
quency dependent) quantity c(co) = — 1/m + G 3 (co). This 
leads to a strong sensitivity of the superconducting DOS 
and T on both u and the chemical potential. Resonant scat- 
tering is usually observed for \c\ <C 1 (corresponding to 
\u\ « °° for constant DOS). In the present case, this trans- 
lates to the condition c(co = 0) <€. 1, allowing for resonant 
scattering even if \u\ <C °°. 

While the case of resonant scattering is of particular 
importance for quantities which are very sensitive to N , 
such as the temperature dependence of the NMR relaxation 
rate, or the London penetration depth, for ARPES exper- 
iments it is not so crucial. The difference between reso- 
nant and non-resonant scattering is most dramatic at small 
impurity concentrations. However, then the self-energy is 
also not so large, and the finite resolution in the ARPES ex- 
periments smears out the sharp feature at resonance which 
exists in the unaveraged spectral function. In order to see a 
clear effect in the EDC, it is necessary to go to rather large 
impurity concentrations n. > 0.02. Hence, in this article, 
we shall concentrate on this case, and consider moderate 
to strong impurity potential strengths u. 

Fig. 3 shows a comparison of the unaveraged spec- 
tral function on the Fermi surface for n. = 0.05 between 
the normal state and the superconducting state at (]) = 
0,45 for various potential strengths. The normal state 
Fermi surface was determined by the criterion that n(k F ) = 
/ oo c/coA(^ f ,co) = 1 /2, i.e., we took the shift of the chemi- 
cal potential caused by the impurity scattering into account 
(of course this shift is only small, of the order of 
for n ( . = 0.05, and dependent on the magnitude and sign of 
u). This criterion ensures that the peak of the normal state 
spectral function is at CO = 0. Concentrating first on the 
spectra for (]) = 45 (at the node), one notes that there is a 
strong dependence of A(k,co) on u. For most values there 
is a clear double peak structure, one peak above, one peak 
below e F . The relative strength of the peaks varies with u. 
As 1 ju is decreased from positive to negative values, the 
peak above E f loses, and the one below e F gains spectral 
weight. At 1/m « 0.0 the two peaks have equal strength, 
and the spectra looks like a single very broad peak. 

This behavior can be readily understood by examining 
Fig. 4 where we plot Im{E n (co)}, which is responsible for 
the quasiparticle damping. We observe a single peak struc- 
ture on a slowly varying background, which moves from 
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a Gaussian broadening of the peaks. Of course this broad- 
ening does not depend on the Fermi surface angle. This is 
very different though for the broadening caused by the fi- 
nite k resolution, as shown in row three. The shape of the 
spectra is similar to a BCS DOS for an isotropic gap, which 
is not surprising, since the order parameter does not vary 
too much within the resolution k window over which one 
integrates. The broadening is clearly strongly dependent 
on k. It is largest at = 45 (the region where the d x2 _ y 2 
gap has a node), and hardly noticeable at = where the 
d x 2_ y2 gap is maximal. This effect is a simple consequence 
of the non-trivial dispersion: Around = 45, the quasi- 
particle band is very dispersive, hence the finite momen- 
tum window of the analyzer samples k points with a wide 
spread in energy. This results in a large broadening. The 
situation is reversed around = ([71,0]), where due the 
saddle point, the band is only weakly dispersive, and hence 
the k points in the resolution window all have similar ener- 
gies leading to extremely small broadening. Note, that the 
k dependence of the OP itself further amplifies this differ- 
ence in momentum broadening, since it has a linear depen- 
dence on k around = 45 (large variation of A fc within 8k), 
and a quadratic one near = (small variation of A k within 
8k). 

Looking at the EDCs calculated according to (1), and 
displayed in the bottom row of Fig. 1, one clearly observes 
that the width of the peaks is about twice as large at = 45 
as compared to = due the momentum resolution effect. 
From this, it is also evident that the widths are k resolu- 
tion limited around = 45, whereas at = they are fre- 
quency resolution limited. The absolute values of the peak 
width are in good agreement with the experimental spectra 
of Ref. 3. This justifies a posteriori the assumption of in- 
terpreting the ARPES data using a simple spectral function 
approach, and also shows that the values claimed by the ex- 
perimentalists about their resolution are self-consistent. 

III. SUPERCONDUCTING GAP 

Next, we are interested in the values of the superconduct- 
ing gap which is obtained by extraction from the calcu- 
lated EDCs. In the experimental literature three different 
methods have been employed: (i) Measuring the shift of 
the leading edge of the EDC in the superconducting state 
as compared to the normal state EDC at the same k point. 
This method is not really suitable to obtain a quantitative 
measure of the gap since the two compared spectra are nei- 
ther taken at the same temperature, nor time. Furthermore, 
the normal state EDCs are substantially broader than the 
superconducting ones leading to an additional uncertainty. 
Therefore, this method can only give an indication that a 
superconducting gap has opened up, but is not trustable for 
an estimate of its size. 

(ii) Measuring the shift of the leading edge of the EDC 
in the superconducting state as compared to an angle- 
averaged spectrum from a non-superconducting reference 



sample such as Ag or Pt which is in electrical contact with 
the superconductor. The spectra of the reference sample 
essentially represent the occupied density of single particle 
states, and should therefore be well represented by a reso- 
lution broadened Fermi function at the temperature of the 
experiment. This method allows for a quantitative estimate 
of the gap, and according to our calculations should always 
give a lower bound of its true value. It is not biased, i.e., no 
assumption about the form of the spectral function is nec- 
essary. However, we shall show that the error as compared 
to the true OP can depend quite strongly on the shape of 
A(k,co). 

(iii) Finally, one can try a fit of the experimental EDC 
to a theoretical curve based on Eq. (1). However, this re- 
quires the knowledge of A(k,co), which is usually assumed 
to be of the simple BCS form. Hence it can only be reliable 
if lifetime effects which could alter the BCS A(k,co) are 
negligible. Furthermore, one has to make an assumption 
about the frequency dependence of the background in the 
spectra which leads to additional uncertainty, even though 
one can check the consistency of the assumptions. Never- 
theless, this is by far the most accurate way to determine 
the gap. 

In this article, we shall use method (ii) and extract the 
shift A sh of our calculated EDCs compared to a resolution 
broadened Fermi function given by N R ((£>,T) = fd(£>'R((£>- 
Co')/((o', T) which is representative for the spectrum of the 
reference sample. The shift is evaluated at half the peak 
value after normalizing the EDC and (co, T) such that the 
maximum of the EDC coincides with the constant value of 
N R ((i>,T) at large negative co. This defines the meaning of 
a normalized spectrum for this article. In the bottom panel 
of Fig. 1, we have illustrated this procedure. Note, that due 
to the finite momentum resolution it is actually possible to 
obtain negative shifts by this procedure, especially if the 
resolution k window is rather small. This is a consequence 
of the non-trivial co dependence resulting from the k aver- 
aging as shown in row three of Fig. 1 . 
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Fig. 2. The superconducting gap on the Fermi surface as a 
function of the angle for the clean case. The true OP A is 
compared with A sh , which is the shift of the leading edge of 
the EDC compared to the reference spectrum iV R (co, T). 

In Fig. 2, we compare the extracted shift with the true 
OP as a function of angle on the Fermi surface. In the 
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a comparison with experiments, and a critical discussion 
will be presented in section V. 

H. SPECTRAL FUNCTION INTERPRETATION OF 
ARPES 

In the analysis of ARPES spectral line shapes, we assume 
the validity of the spectral function interpretation. In other 
words, we assume that the ARPES intensity (neglecting ef- 
fects of finite experimental resolution) is given by 7(k, co) = 
7 (k)/(co)A(k,co). The pref actor I Q (k) depends on incident 
photon energy, polarization, and on the final state through 
the electron-photon matrix element, however, it is only 
weakly dependent on CO or temperature. Hence, the line 
shape is entirely determined by the Fermi function /(co), 
and A(k,co). We shall later comment on this simplifying 
assumption, and show that it leads to a self-consistent in- 
terpretation of the ARPES experiments on Bi 2 Sr 2 CaCu 2 8 
(Bi2212), which is the best studied cuprate. Recently, it 
has been demonstrated that this approach works also well 
in the analysis of other quantities. 9 

A slight complication regarding the interpretation of 
experimental spectra results from the finite resolution of 
the detector in both frequency and momentum. Usually, 
the detector probes a finite momentum window 8k with 
constant probability distribution due to its finite extent. 
Apart from the physical dimension of the detector, the size 
of this window depends also on the incident photon energy. 
A circular window of radius 0.036ft/ a, a the planar lat- 
tice constant, is a typical dimension, realistic for the mea- 
surements from Ref. 3 performed at incident photon en- 
ergy 22eV. All calculations in this paper were done using 
this k resolution, and the temperature entering the Fermi 
function was fixed at T = IQK, a typical value for actual 
measurements in the superconducting state. The frequency 
resolution of the detector is a Gaussian, and current state- 
of-the-art devices can reach a standard deviation a as low 
as 7meV. Taking both resolution effects into account, the 
ARPES intensity or EDCs should be given by 3 

7(k,co) =I J^dk' J dco'fl(co-co')/(co')A(k',co'), (1) 

where R((ti — co') ~ exp ^— ^ m 2 ~ a m - j . Here we neglected the 
k dependence of I Q . 

It is very instructive to calculate EDCs according to this 
formula within simple BCS theory. Then the spectral func- 
tion takes the familiar form 

A(k,co) = M 28(co-£ k ) + v2S(co + £ k ), (2 ) 

with E k = y^~+A£ the quasiparticle energy, and w k = ( 1 + 

£ k /E k )/2, v k = (1 - ^jE k )/2 the coherence factors. To 
be as realistic as possible, for the quasiparticle dispersion 
^ k , we use a tight-binding fit to normal state ARPES data 
on Bi2212 3 ' 10 with real space hopping matrix elements 



0' 



[0.879,-1,0.28,-0.087,0.094,0.087], (t on- 



site, /j nn, t 2 nnn hopping, ... ). n ' 12 All energies are mea- 
sured in units of \t l \ = 0.149eV. The value of t Q corresponds 
to a hole doping of 8 = 0.17. As suggested by the ARPES 
data, the fitted quasiparticle band-structure exhibits a sad- 
dle point at k = (71,0) with the resulting van Hove singular- 
ity located at approximately 30meV below e F in Bi2212. 
Normal state data on other hole-doped cuprates show a 
similar dispersion 
1.0 
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Fig. 1. The effect of the finite experimental resolution on the 
EDCs as calculated from the pure BCS d x2 _ y2 spectral func- 
tion. The three columns are calculations for different angles 
on the Fermi surface, ij) = 45, 35,0 from left to right. Top row: 
the pure spectral function. Second row: spectral function after 
frequency broadening. Third row: spectral function after mo- 
mentum broadening. Last row: EDC calculated from (1), and 
normalized to the spectrum of a resolution broadened Fermi 
function (dotted curve). 

The OP is chosen as A k = A(cosfe^ — cosfc^) /2, and the 
size of the maximum gap is fixed at A = 0.2 (for T = 
0), which corresponds to approximately 30meV, consistent 
with ARPES results on Bi2212. ' Fig. 1 demonstrates 
the effect of the experimental resolution on the ARPES 
line shape. We show spectra at three different angles (]) = 
0, 35 , 45 on the Fermi surface. (]) is measured relative to the 
(jc,7t) point, around which the Fermi surface is closed. The 
top row displays the pure BCS spectral function, i.e., two 
delta functions provided the OP is finite (the node in the 
d x 2_ y 2 gap is at (j) = 45). In the second row, the effect of 
the frequency resolution is shown, which simply leads to 
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Effect of non-magnetic impurities on the gap of a d x2 _ y2 superconductor 
as seen by angle-resolved photoemission 

R. Fehrenbacher 
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An analysis of angle-resolved photoemission (ARPES) experiments in the superconducting state of the high 
T copper-oxides is presented. It is based on a phenomenological weak-coupling BCS model which incorpo- 
rates the experimental normal state dispersion extracted from ARPES, and non-magnetic impurity scattering 
in the presence of a d x2yl order parameter (OP). It is shown that already in the pure case, the broadening 
by finite momentum resolution of the analyzer leads to a finite region of apparent 'gaplessness' around the 
true node of the OP. Non-magnetic impurities further amplify this effect by introducing additional spectral 
weight around zero frequency. At sufficiently large impurity concentrations n. sa 0.02 — 0.05, this results in 
an extended region of 'gaplessness' up to S(|) = ±7 (0 the angle on the Fermi surface) around the true node 
for a large range of moderate to strong impurity potential strengths. Different ways to identify the presence 
of impurity scattering in the ARPES spectra are proposed. 
PACS numbers: 74.20.-z 74.20.Mn, 74.25.Nf, 74.62.Dh 



I. INTRODUCTION 

Major experimental progress has recently been made in 
identifying the order parameter (OP) symmetry of the 
high T cuprates. A larger number of results strongly 
suggest that the carriers pair with a non-trivial, spatially 
anisotropics? 2 _ 2 symmetry. 1 Among these probes, angle- 

x y ^ 

resolved photoemission (ARPES) has the unique potential 
of determining the superconducting gap as a function of 
momentum, and indeed, strong support for the d- wave sce- 
nario has recently been reported using this technique. 2 ' 3 

Since the cuprates are the first materials in which 
ARPES experiments were able to detect a superconduct- 
ing gap, there is not much theoretical work experimental- 
ists can compare their spectra to, and hence in most cases 
the extraction of the gap is done on a rather qualitative 
level. It is therefore important to provide some minimal 
theoretical guidance, and put the data analysis on firmer 
grounds. To this end we calculate so-called energy distri- 
bution curves (EDC) using a rather simple phenomenolog- 
ical BCS model which consists of the ARPES normal state 
quasiparticle dispersion supplemented by a d 2 _ 2 pair in- 

x y 

teraction and non-magnetic impurity scattering. 4 The ef- 
fects of impurity scattering are of particular importance, 
since the reduced symmetry of the OP allows for severe 
pair-breaking even in the non-magnetic case. 5 

One rather selective signature of the of-wave state, as 
compared to a strongly anisotropic s-wave state for in- 
stance, is the possibility of so-called resonant scattering 
by non-magnetic impurities. 6 ' 7 If the impurity potential is 
such that the resonant condition is satisfied, a very small 
impurity concentration n. is sufficient to create a large 
residual density of single particle states (DOS) A^ res at the 
Fermi level. Naturally, this effect also has to manifest itself 
in the spectral function A (k,co), and hence it is of interest 



to study whether it could be observed in an ARPES exper- 
iment. Even away from resonance, the impurities induce a 
finite N , albeit the growth with the number of impurities 
is much smaller. Nevertheless, at sufficiently large concen- 
trations of a few percent, one expects to see an effect on 
A(k, co) in this limit also. In earlier work, 8 we already stud- 
ied the influence of impurity scattering on ARPES spectra 
from a slightly different perspective, but the model used 
there did not incorporate a realistic normal state quasipar- 
ticle dispersion. When comparing to experimental spec- 
tra, however, a realistic model for the dispersion is crucial, 
since it dramatically affects resolution broadening, as we 
shall discuss below. 

By analogy with the rise of N as a function of impu- 
rity concentration, one might expect that the effect on the 
gap as extracted by ARPES is also to progressively wipe it 
out as impurities are added. Indeed, we shall show that, de- 
pending on how the gap is extracted from the ARPES spec- 
tra, there can be a finite region on the Fermi surface around 
the d x 2_ y2 node, where the measured gap appears to van- 
ish, even though the order parameter is finite. However, to 
a large extent, this effect is caused by the finite momentum 
resolution of the analyzer. It occurs even for the case of a 
pure BCS spectral function, although the 'gapless' region 
is then much smaller. 

The plan of this article is as follows: In section II, we 
shall start with a brief discussion of the photoemission pro- 
cess, and discuss the influence of finite experimental reso- 
lution in both energy and momentum on the spectral line 
shape. Based on this, in section III, we shall discuss dif- 
ferent methods to extract the superconducting gap from the 
ARPES spectra, and show how well its resulting k depen- 
dence compares with the true OP. In section IV, the influ- 
ence of non-magnetic impurity scattering on the line shape, 
as well as on the extracted gap will be elucidated. Finally, 
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